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ABSTRACT 



Aims. We investigate the spot activity of the young chromospherically active main sequence star LQ Hya. Our aims are to identify 
possible active longitudes, estimate the differential rotation and study long and short term changes in the activity. 
Methods. Our analysis is based on 24 years of Johnson V-band photometry of LQ Hya obtained with the T3 0.4m Automated 
Telescope at the Fairbom Observatory. We use the previously published Continuous Period Search (CPS) method to model the evolu- 
tion of the light curve of LQ Hya. The CPS fits a Fourier series model to short overlapping subsets of data. This enables us to monitor 
the evolution of the light curve and thus the spot configuration of the star with a higher time resolution. 

Results. We find seasonal variability in the mean level and amplitude of the light curve of LQ Hya. The variability of the light curve 
amplitude seems not to be cyclic, but the long-term variations in the mean magnitude could be explained by an approximately 13 year 
cycle. However, because of the limited length of the observed time series, it is not yet possible to deteiTnine whether this structure 
really is periodic and represents an activity cycle. Based on fluctuations of the light curve period, we estimate the differential rotation 
of the star to be small, and the star is potentially very close to a rigid rotator We search for active longitudes from the inferred epochs 
of the light curve minima. We find that on time scales up to six months there are typically one or two relatively stable active areas 
on the star with limited phase migration. On the other hand, on time scales longer than one year, no stable active longitudes have 
been present except for the period between 2003 and 2009 and possibly also some time before 1995. Neither do we find any signs of 
flip-flops with a regular period. The mean time scale of change of the light curve during the observation period is determined to be of 
the same order of magnitude as the predicted convective turnover time for the star. 

Key words. Stars: activity, starspots, individual: LQ Hya 



1. Introduction 

LQ Hya (HD 82558) is a young single chromospherically active 
st ar (V ^ 7.8 B- V = 0.9, K2V) classified as a BY Dra star 
by Fekel et alJ(ll986bl) . Its strong activity is clearly evident from 
the su bstantial Ca HK emission, logTJ^j^ - -4.06 JWhite et alj 
l2007h . p lacing it confidently within the "very active" regime de- 
fined bv lHenrv et aP (Il996h. 

It was suggested bv lFekel et al.l ( Il986al) that the star is a very 
young obj ect which has recen tly arrived at the zero-age main- 
sequence. iMontes et al.l (l200lh classified it as a member of the 
young disc popu lation following the d e finitio n of lEggenI (Il984bl 
119891) . Recently 'Nakai ima & Morinol (12012*) identified the star 
as a possible member of t he IC 2391 superclus ter thus estimating 
its age to be 35-55 Myrs dMontes et al.ll2001h . 

The magnetic activity of LQ Hya is strongly manifested as 
starspo ts causing rotational modulation of brightness (Eg gen 
ll984 aV The rotation peiiod of the star can be measured from 
this modulation. Previously the rotation period of LQ Hya has 
been estimated as P = 1.601 136 ± 0.000013 d by Ij etsu ( 1993|), 
P = 1.601052 + 0.000014 d by 'Berdvugina et aL (2002) and 
P = 1.60066 + 0.00013 d by Kovari et al. (2004). Some vari- 
ability of the observable photometric rotation period is expected. 



* The analysed photometry and numerical results of the 

analysis are both published electronically at the CDS via 

anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A-HA/yyy/AxxX| 



however. This can be caused by the spots moving with different 
angular velocities governed by the underlying surface differen- 
tial rotation or large scale magnetic field. Also changes in the 
light curve shape governed by active regions growing and decay- 
ing at different locations on the stellar surface can cause period 
variations unrelated to the spot movements. 

The diff erential rotat ion of LQ Hya is found by many authors 
to be small. iJetsul (Il993h analysed the fluctuations of the photo- 
metric period of LQ Hya within 3cr limits and reported the rela- 
tive scale of them to be Z w 0.015. This value can be interpreted 
as the relative scale of rotation periods of the observed starspots. 
Provided that the spots trace the surface rotation of the star and 
that there have been spots on all stellar latitudes from the equator 
to the poles, it thus gives an estimate for the differential rotation 
coefficient k, i.e. indicating k ^ 0.015. You (2007) used a similar 
approach to estimate the differential rotation of the star from the 
light curve period fluctuations and derived the value k x; 0.025. 
Both of these values are also similar to the theoretical estimate 
k ^ () .0128 obtained for LQ Hya bv iKitchatinov & Olemskovl 
(1201 ih using a mean field hydrodynamical model. 

Even smaller differentia l rotation values we re reported by 
iBerdvugJna et al.l ([2002) and lKovari et al.l (|2004 who both used 
Doppler image s of LQ Hya in c o njunc tion with photometry in 
their analyses. 'Berdvu gina et al.l (l2002l) compared photometric 
periods of LQ Hya to Doppler images coiTesponding to the 
same epoch of time to infer the latitude of the main spot. Using 
this approach they retrieved a differential rotation coefficient of 
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k K 0.002. iKovari et alJ (l2004l) used cross correlation between 
adjacent Doppler images to estimate the differential rotation and 
reported A: = 0.0057. 

Many authors have reported sign s of cyclic behaviour in the 
spot activity of LQ Hva. lJetsul(ll993h found a 6.24 yr cycle in the 
mean brightness of the sta r using light curve fits for 9 years of 
photometric observations. ' Olah et aP (1200 0*) applied a Fourier 
analysis to 16 years of photometry and found cycle periods of 
both 6.8 yr and 11.4 yr in the mean brightnes s. Based on light 
curve inversions from 20 years of photometry, Berdvugina et al.l 
(I2OO2I) reported a 15 yr cycle for the mean brightness, a 7.7 yr 
cycle for the light curve amplitude modulation and a 5.2 yr flip- 
flop cycle. They a lso identified a 1:2: 3 resonance between the 
cycle frequencies. iKovari et al.l (|2004|) used Fourier analysis to 
search for cycles in the light curve mean brightness from 8 years 
of photometry and found a possible cycle of 13.8 yr along with 
its first harmonic 6 .9 yr and also wea k signs for a 3 .7 yr cycle. 

In addition to iBerdyugina et al.l (.2002.) . also iJetsul (il993h 
reported active longitudes for LQ Hya. The nature of 
the active longitudes was, however, quite different. Where 
[Berdvugina et al. (2003) r eported tw o active longitudes with 
A<p - 0.5 phase separation. iJetsul (1 1 99 3 ) claimed the phase sep- 
aration between the active longitudes to be only A(p - 0.25. 

Doppler imag es of L Q Hya have been r econstru cted by 
[Strassmeieretal. ('1993'), Rice & Strassmeier' ('199J||, iDonatil 
( 1999), Donati et al. (2003) and Kovari et al. (2004). They typi- 
cally show spotted areas at mid latitudes relatively far from the 
visible pole. From time to time, there have been longitudinal 
concentrations of spots in these maps, but no clear pattern of sta- 
ble active longitudes. Occasionally the surface reconstructions 
have shown a complete latitudinal band of spots encircling the 
star On the other hand, reconstructio ns of the surface magnetic 
field u sing Zeeman Doppler imaging ( IDonatill 19991: iDonati et al.l 
I2OO3I) have sometimes shown opposite magnetic polarities on 
different sides of the visible pole. This may indicate at least an 
occasional presence of active longitude like features on the star 



2. Analysis of the data 

The photometry of LQ Hya was obtained over a 24 year 
time span between HJD - 2447141 (11 December 1987) and 
HJD = 2455684 (3 May 2011) with the T3 0.4 m Automatic 
Photoelectric Telescope (APT) at the Fairbom Observatory in 
Arizona. The complete dataset analysed in this paper consists of 
267 1 Johnson V-band observations of LQ Hya minus the com- 
parison star HD 82477 (V = 6.13, B-V ^ 1.19). To monitor 
the constancy of the comparison star, 2272 additional simultane- 
ous V-band observations of the check star HD 82428 (V = 6.14, 
B -V - 0.24) were obtained. 

The typical error of the target star photometry is estimated 
to be between 0.003 a nd 0.004 magnitudes from observations of 
constant stars (lHenrvl fl995). The errors of the check star obser- 
vations are expected to be somewhat larger as fewer individual 
integrations have been used to determine their values. For a brief 
des cription of the operatio n of the APT and reduction of the data, 
see lFekel & Henrvr(l2005h and references therein. 

Our Johnson V band LQ Hya minus comp star and check star 
minus comp star differential magnitudes are presented in the top 
and bottom panels of Fig. [T] respectively. LQ Hya shows signifi- 
cant variability on night-to-night (rotational), year-to-year (spot 
lifetime), and decadal (magnetic cycle) time scales. The check 
star observations, plotted on the same scale as the LQ Hya ob- 
servations show no evidence for variability on any time scale 
in either the comp star or the check star The somewhat larger 



1.65 




. .:; 




;!; 












;. P 




I 


1 l» 




"i ' " 


' ^ ■ '■: 


:■>■■ ; 


■ i^ 






■ 


■■ ■ *■■! 


V. 1^ 


-'■ 


■J: :V 


1.75 


y 


f ' : 


■ '■''■ i- 
:''| 


%} 


m 


i^ 


• 




i\ » i- 

;■■ -i :: 

:■. V ;;; 




y 


|| 


1.80 


'■' iS' 














■., ?' ■' 






. 




■ ;'-' 












■■■■i 


' : .-V ■■■: 


■ .:'= 










;*; 
























1.85 


, «■ 














\ ■ 


"' ■. 


:^ 




" 


1.90 




















■. 




- 



0.00 

0.02 
0.04- 



1- 



n«?si^^ll^sl^^lsiii>Sii 



2000 

year 



Fig. 1. Top: V band differential photometry of LQ Hya minus 
comparison star The small error bar in the upper right comer 
denotes the scale of the photometric uncertainty of ±0.004 mag. 
Bottom: V band differential photometry of the check star minus 
comparison star. The scale is the same in both of the panels. 



scatter and slight offsets seen in the first few years of the check 
minus comp star observations are the result of instrumental up- 
grades that affect mainly the check minus comp differential mag- 
nitudes because of the large color difference (A(B -V) - -0.95) 
compared to LQ Hya - comp (A{B -V) - -0.29). 

The photometry of LQ Hya was analysed using 
the Continuous Per iod Search (hereafter CPS) method 
(iLehtinen et al.l I2OIII) . The method models the photometry 
in short subsets of the data using a low Kth order Fourier series. 
This modelling provides estimates for the mean differential 
magnitude M, peak-to-peak amplitude A and period P of the 
light curve, as well as epochs for the primary and a possible 
secondary minimum fn,in,i and f,nin,2- To allow for variability 
in the model complexity, the CPS performs fits using Fourier 
series of orders K = . . .2. The model order K = Q corresponds 
to a simple constant brightness model and describes the absence 
of any intrinsic variability in the light curve. For each dataset, 
the Bayesian inform ation criterion is app lied to determine the 
best modelling order (ILehtinen et al.ll201ll Eq. 6). 

The period search was done within a ±5% interval around 
the a priori period estimate Pq. The search was limited to this 
interval because of the risk for interference with spurious peri- 
ods. In line with the previous results from LQ Hya we chose the 
value of Po = 1 -6 d. The upper limit for length of the individual 
analysed datasets was set at Ar^ax = 30.4 d, which is 19 times 
the length of Pq. To get good time resolution for the evolution 
of the light curve parameters [M,A, P, fmin.i, fmin,2], the datasets 
overlap with each other so that a new dataset was selected af- 
ter 1 d from the start of the last one. Only datasets containing 
n > 12 observations were included in the analysis. The mean of 
the residuals of all the model fits is ? = 0.009. 



3. Results 

The CPS automatically divides the data into segments consist- 
ing of mutually overlapping datasets. In addition to this, it cleans 
the data by removing outl iers and temporally isolated data points 
(see Le htinen et al.ll20I 11) . The segment division for the data of 
LQ Hya is summarised in Table [1] The table lists the dates of 
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Fig. 2. Light curves and light curve fits of the 111 independent datasets. The ephemerides for the phases of the observations and 
light curves are explained in Sect, |3] 



the first and last analysed data point in each segment, the lengths 
of the segments in days, as well as the numbers of all analysed 
datasets and independent datasets in the segments. The provi- 
sional segment SEG 5 is lacking from the listing because it did 
not contain enough data for modelling, i,e, it contained only 
some isolated data points. Overall, the data from the first few 
years is fragmentary, resulting in some shorter segments. 

The mutually independent datasets were selected from all 
datasets so that they do not overlap with each other Using the 
complete set of analysed datasets gives a detailed view into the 
time evolution of the light curve. On the other hand, using only 
the independent datasets removes any effects introduced by mu- 
tual correlations between the models of partially overlapping 
datasets. 

The numerical results for all of the 1077 analysed datasets 
can be accessed electronically at the CDS, Light curve fits for the 
111 independent datasets are presented in Fig, |2] Observations 
in the datasets having a minimum epoch fmin,! and period P have 
been folded to phases (p - FRAC[(f - fmin.i)/^], where FRAC[jt:] 
removes the integer part of x. To track the migration pattern of 
the primary minima, the phases have been shifted with the phase 
of the primary minimum 0ai computed using the ephemeris of 
Eq, [3] and an additional phase shift of (pshift - 0,5, to make the 
minima more apparent. The ephemeris of Eq, |3]is defined later 
in Sect, l3,3l in the context of active longitudes. 



3. 1 . Long term variability of M, A and P 

The long term evolution of the light curve mean M, amplitude 
A and period P are presented in Fig, |3] We find both regular 
trends and random fluctuations from these light curve parameters 
during the last 24 years of observations. 

The regular variations are most striking in the light curve 
mean M which naturally follows the overall variation of the raw 
photometry (Fig, [T]i, The trend in the M variations resembles 
what may turn out to be a regular activity cycle. There was a 
brightness minimum at 1989 and again around 2002, Between 
these two minima the mean brightness of the star increased with 
an amplitude of nearly 0,1 mag. After the 2002 minimum the 
mean brightness has climbed steadily towards a new maximum. 

The variations in M could be explained with an activity cycle 
of approximately 13 years. However, the total length of observa- 
tions only includes what seems to be one and a half cycles, so 
conclusively proving the existence of such a cycle, let alone ac- 
curately determining its length, remains unfeasible. Such consid- 
erations are rendered even harder due to the fact that the profiles 
of these two apparently separate cycles in M are quite different. 
What may appear as cyclic behaviour in the light of the available 
data may turn out to be random variability in the future and vice 
versa. 

The light curve amplitude A has varied between 0,00 and 
0,24 mag. The short-term variations are quite regular, but do not 
show any suggestive signs of cyclic behaviour In particular, cor- 
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Table 1. Summary of the segments: Dates of the first and last 
analysed data point in HJD and calendar dates (yy/mm/dd), seg- 
ment length ArsEG rounded into full days, the number of indi- 
vidual datasets nsET and the number of mutually independent 
datasets «ind- 



SEG HJD -2400000 date 



ArsEG "SET «IND 



1 


47199- 


- 47230 


88/02/08 - 


- 88/03/10 


31 


2 


1 


2 


47277 - 


- 47304 


88/04/26 - 


- 88/05/23 


27 


1 





3 


47460 - 


- 47660 


88/10/25 - 


- 89/05/14 


201 


69 


6 


4 


47832 - 


~ 48027 


89/11/02- 


- 90/05/16 


195 


47 


5 


6 


48348 - 


- 48394 


91/04/02- 


-91/05/18 


46 


7 


1 


7 


48696 - 


- 48759 


92/03/07 - 


- 92/05/17 


63 


6 


1 


8 


48911- 


-49132 


92/10/15 - 


- 93/05/25 


222 


32 


4 


9 


49290 - 


_ 49499 


93/10/30 - 


- 94/05/27 


209 


60 


6 


10 


49645 - 


~ 49866 


94/10/19 - 


- 95/05/29 


222 


49 


6 


11 


50006 - 


- 50226 


95/10/15 - 


- 96/05/23 


221 


67 


5 


12 


50391 - 


- 50582 


96/11/04- 


- 97/05/14 


191 


54 


4 


13 


50736 - 


- 50955 


97/10/14 - 


- 98/05/22 


220 


61 


6 


14 


51103- 


-51325 


98/10/16 - 


- 99/05/27 


223 


93 


7 


15 


51474- 


-51687 


99/10/23 - 


- 00/05/23 


213 


69 


6 


16 


51861- 


- 52052 


00/11/13- 


-01/05/23 


191 


20 


5 


17 


52287 - 


-52421 


02/01/12 - 


- 02/05/27 


135 


36 


4 


18 


52594 - 


- 52785 


02/11/15- 


- 03/05/23 


192 


46 


5 


19 


52977 - 


-53149 


03/12/04 - 


- 04/05/24 


172 


37 


5 


20 


53329 - 


-53506 


04/11/20- 


- 05/05/16 


177 


37 


4 


21 


53660 - 


-53876 


05/10/16 - 


-06/05/21 


217 


43 


5 


22 


54044 - 


- 54238 


06/11/05- 


-07/05/18 


194 


55 


6 


23 


54400 - 


- 54599 


07/10/26 - 


-08/05/13 


200 


57 


6 


24 


54810- 


- 54966 


08/12/10 - 


- 09/05/15 


156 


26 


4 


25 


55126- 


-55301 


09/10/21 - 


- 10/04/15 


176 


31 


4 


26 


55499 - 


- 55684 


10/10/29 - 


-11/05/03 


186 


72 


5 



relation between the independent M and A estimates is absent 
with the linear Pearson correlation coefficient being r - 0.10. 

When the amplitude is at its lowest, the data is best described 
by a constant brightness model. This indicates that occasionally 
the light curve of LQ Hya reduces to such low amplitudes that 
the periodic variability is buried under the random observational 
errors and systematic errors introduced by the model. Physically 
this means a nearly total absence of spots or alternatively an ax- 
isymmetric spot distribution. On the other hand, the maximum 
values of A above 0.2 mag are a substantial indication of strongly 
concentrated spot activity. 

The variability of the light curve period P seems to consist 
of random fluctuations. These could be caused by starspots ran- 
domly occurring on diff'erent latitudes having different rotation 
periods and thus tracing differential rotation. Alternatively they 
may be caused by changes in the light curve shape as some active 
regions decay and others form at different longitudes. Following 
the assumption that the random fluctuations are caused by dif- 
ferential rotation, they are used to give an estimate for the differ- 
ential rotation of LQ Hya in Sect. 13.21 

3.2. Differential rotation 

An estimate for the differential rotation of LQ Hya was obtained 
from period fluctuations using the formula (|Jetsulll993h 



Z = 



6AP„ 



(1) 



where w,- - CTp^ are the weights, -Pw = H w/f ,/ 2 >v/ is 
the weighted mean and AP^ = VS '<'^S.Pi - ^w)^/ Z w'/ is the 
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Fig. 3. Long term variation of the mean differential magnitude 
M (top), total light curve amplitude A (middle) and photometric 
period P (bottom). M and A are given in magnitudes and P in 
days. Squares with error bars indicate parameter estimates from 
the independent datasets and dots parameter estimates from all 
other reUable datasets. 



weighted standard deviation. The parameter Z measures the rel- 
ative variability of the light curve period within its +3cr limits. 

When estimating differential rotation from photometric pe- 
riod fluctuations, the range of period fluctuations is often equated 
directly with the differential rotation coefficient k. This can gen- 
erally not be assumed. The period fluctuations can provide infor- 
mation of the rotation period only from that latitude range of the 
stellar surface on which the spot activity occurs. As this range 
might be quite limited, we do not expect the photospheric period 
fluctuations to correspond to the total range of photospheric ro- 
tation periods present on the star. However, in the absence of any 
knowledge about the true latitude range of the spot activity, we 
may expect k to be of the order of Z or somewhat greater 

For LQ Hya we get P^ ± AP„ = 1.6043 + 0.0052 d, coiTe- 
sponding to fc « Z = 0.020 or AQ w 0.078 rad /d. This is in line 
with previous estimates, especially those by iJetsul (Il993h and 
You (.2007.) . who estimated the differential rotation in a similar 
manner. 

We add a few caveats to the interpretation of the pe- 
riod fluctuations as tracers of differential rotation. First, if the 
starspots are caused by a large scale dynamo field, it is pos- 
sible that they do not fol low the surface rotation of the star 
jKorhonen & Elstnenl201 lb . It is also possible that the small pe- 
riod fluctuations are not caused by starspots having different ro- 
tation periods around the star but rather due to starspot growth 
and decay affecting the light curve shape. Lastly, we note that in 
datasets that have few data points and a low light curve ampli- 
tude there is a certain amount of ins tabiUty in the period detec- 
tion (Table 2 in lLehtinen et al.ll2011h . 

3.3. Minimum pliases ^min snd active longitudes 

The longitudinal distribution of the spots on LQ Hya can be stud- 
ied using the epochs for the primary and secondary minima t^i„i 
and fmin,2 of the modelled light curve. For any reasonable rota- 
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Fig. 4. Primary (black) and secondary (grey) light curve minimum phases, 0min,i and <^min,2, for segments that have Wiei > 10 datasets 
with reliable light curve models. The minimum phases are computed from the primary and secondary light curve minimum epochs, 
fmin.i and fmm,2, with the median light curve periods Pmed of each segment and the first primary minimum set at = 0.25. None of 
the secondary minima are considered fully reliable in the sense that they appear in less than 95% of the bootstrap runs performed 
by the CPS. 



tion period estimate P, these can be transformed into light curve 
minimum phases 



= FRAC 



(2) 



where FRAC[x] removes the integer part of x. 

The minimum phases ^min,! and 0niin,2 for individual seg- 
ments with «i-ei > 10 datasets with reliable parameter estimates 
are presented in Fig. |4] The phases are folded from the primary 
and secondary minimum epochs using the median photometric 
period Pmed of cach segment and placing the first primary min- 
imum of each segment at = 0.25. For each segment, time is 
given starting from the the first analysed data point in the seg- 
ment as given in Table [T] 

A recurring pattern between the segments has been the pres- 
ence of one or two active regions wandering slightly in phase 
(Fig.|4|. In many segments the two active regions inferred from 
the light curve minima stayed roughly at the opposite sides of 
the star resembling long-lived active longitudes. In other seg- 
ments, such as SEG 14 and SEG 26, there are clear examples 
of one light curve minimum splitting into two or two minima 
merging into one as the underlying active regions have moved 
away or towards each other At some critical phase separation 
the two active areas have moved too close to each other to pro- 
duce separate light cur ve minima and instea d we only observe 
one merged minimum (JLehtinen et aLllMTTl Eq. 12). In several 
segments (SEG 3,10, 15, 16 and 26) the main photometric min- 
imum switched from one phase to another Although occurring 
only on a short timescale, this may be analogous to the flip-flop 
behavior discovered in other stars (e.g. Jetsu 1996). 

To examine whether there have been any long lasting ac- 
tive longitudes on LQ Hya, i.e. that the light curve minima have 



clustered arou nd certain phases, we performed the Kuiper test 
(iKuipeJ 19601) for the minimum epochs t^nm- Our formulation for 
the test is from Ijetsu & Pelt (1996) . This procedure computes 
the Kuiper periodogram for tmm and tests the null hypothesis 
of uniform (i.e. random) phase distribution. It also determines 
the critical levels for the most significant periods. Examples of 
apphcation of this m ethod can be found in iJetsul (Il996h and 
lLehtinenet"ai] (120111) . 

The Kuiper test was performed for all reliable light curve 
minimum epochs from the independent (i.e. non-overlapping) 
datasets. This set comprised of 1 1 1 primary minimum epochs. 
Although many of the model fits provided by the CPS contained 
also secondary minima, none of these were considered reliable 
in the sense of being present in more than 95% of the bootstrap 
runs performed by the CPS ( Lehtinen et al.ll2()ITI) . Consequently 
none of the secondary minimum epochs fmin.i were included in 
the Kuiper test. 

The most significant period found for the reliable indepen- 
dent minimum epochs was P^i = 1.603693 + 0.000058 d with 
the critical level 2 = 6.1- 10"^. All independent minima, folded 
into minimum phases (p^in with this period, are presented in 
the top panel of Fig. |5] The structure responsible for this pe- 
riod is clearly visible between the years 2003 and 2009. During 
this time interval the primary light curve minima were confined 
within a tight phase region. There are no signs of this structure 
being extant either before 2003 or after 2009, where the mini- 
mum phases show no structure with this period. An approximate 
ephemeris for the central meridian passing of this active longi- 
tude is given by 



HJDmin = 2447201.3 -i- 1.60369£. 



(3) 
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Even between the years 2003 and 2009 the active longitude has 
not remained stable but has moved back and forth within a phase 
range of A0 a; 0.2. A secondary active longitude, consisting 
solely of secondary light curve minima, appears to have been 
present for some years after 2003 with a phase separation of 
A^ «j 0.5 from the primary active longitude. 

We performed additional Kuiper tests for the light curve min- 
imum epochs to investigate the minimum phase distribution be- 
fore and after the emergence of the obvious active longitude dur- 
ing 2003. The analysis was done for two samples of independent 
reliable minimum epoch estimates. The first sample consisted of 
the 67 independent primary minimum epochs in segments SEG 1 
- SEG 17 (i.e. until April 2002), while the second sample con- 
tained the remaining 44 independent primary minimum epochs 
in segments SEG 18 - SEG 26 (i.e. starting from November 
2002). 

For the latter part of the minimum epochs, the analysis found 
the period -Pai,i8-26 = 1.60374 + 0.00013 with the critical level 
Q = 2.3 • 10"^. This corresponds to the same periodicity than P^i 
found for all minima but now with a significantly lower critical 
level. The minimum phases folded with this period (ephemeris 
HJDmin = 2447201 .3 + 1 .60374£') are shown in the middle panel 
of Fig. |5] and follow very closely the pattern of the minimum 
phases folded with the period P^\. 

For the first part of the minimum epochs the analysis indi- 
cated PaU-17 = 1.68929 + 0.00008 d to be the most signif- 
icant periodicity with the critical level Q - 7.7 ■ 10"^. The 
minimum phases folded with this period (ephemeris HJDmin - 
2447200.8 + 1.68929£') are presented in the bottom panel of 
Fig. |5] This folding brings out some structure for the primary 
minimum phases during the year 1995 and some time before 
that. However, the period Pa\.\-\i is considerably longer than ei- 
ther of the other two active longitude period estimates P^x and 
^ai,i8-26 or even the highest individual light curve period esti- 
mate Pmax = 1.6297 d. It is thus possible that the period Pai,i-i7 
corresponds rather to a spurious effect than a real active longi- 
tude. 

The rest of the minimum epochs seem to exhibit no periodic 
structure at all. 

3.4. Time scale of change 

In addition to determining the light curve parameter estimates, 
the CPS also determines the time scale of change Tq for each 
individual dataset. This is defined as the time during which each 
model fit remains an adequate descripti on for the subseq uent 
datasets. As we have demonstrated (Le htinen et al.l l201ll) . the 
value of Tc can change dramatically from one dataset to another 
A more meaningful value is the mean time scale of change Tc, 
which is just the mean of all Tc. 

For LQ Hya we get Tc = 50.5 d based on the individual Tc 
values of all independent datasets. This is longer than the maxi- 
mum length of the datasets Armax - 30.4 d by a wide margin. In 
other words, the light curve of LQ Hya typically retains its shape 
during the whole dataset and our choice for the dataset length is 
well justified. 

As the value of Tc estimates the typical time scale associ- 
ated with the evolution of the spot configuration, it is interest- 
ing to check how this compares with the theoretical convective 
turnover time t^. . Previously, we used the interpolation formula 
bylOssendriiver ( 1997) 

Tc = -68.3 + 224.8(B - V) - 177.2(B - Vf + 57.0(B - V)\ (4) 
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Fig. 5. Light curve minimum phases folded according to the 
ephemeris HJD™,, = 2447201.3 + 1.60369£ for all of the min- 
imum epochs (top), HJDmi„ = 2447201.3 + 1.60374£' for the 
minimum epochs starting from November 2002 (middle) and 
HJDmin = 2447200.8 -i- 1.68929£' for the minimum epochs until 
April 2002 (bottom). Black squares denote primary minima and 
grey triangles secondary minima. 



to estimate the convectiv e turnover time for the young so- 
lar analogue HD 116956 (" Lehtinen et 10 12011*). This formula 
is based on the theoretical calculations of Kim & Demarqi^ 
( 1996). By using the H ipparcos value B-V = 0.933 for LQ Hya 
(.Ferryman et al.lfl997l) . we get r^ = 33.5 d. Similarly as in the 
case of HD 1 16956, the values of Tc and r^ are of the same order 
of magnitude, t^. being about 65% of the length of Tc. 



4. Conclusions 

We have analysed 24 years of Johnson V-band photometry of 
the chromospherical ly active star LQ Hya with our recently pub- 
lished CPS method ("Lehtine n et al.ll201l1) . This method models 
the observed photometry with a variable order Fourier serie s us- 
ing a sliding window for choosing the analysed datasets. The 
modelling provides estimates for the light curve mean differen- 
tial magnitude M, total light curve amplitude A, photometric pe- 
riod P and primary and secondary light curve minimum epochs 
fmin.i and fmin,2 as functions of time. 

LQ Hya has displayed regular variability in the M and A es- 
timates. Especially the variability of the mean differential mag- 
nitude M resembles a segment from a quasi periodic time se- 
ries. The variability could be explained with a roughly 13 year 
spot cycle. This is undoubtably the same structure that is be- 
hind the 11.4 yr c ycle reported b y Olah et al. (2000), the 15 yr 
cycle reported by Berdyuginaet al. (2002) and the 13.8 yr cy- 
cle suggested by Kovari et al. (2004). However, our present data 
only includes one and a half cycles of this suggested activity cy- 
cle. This is not enough to conclude whether the mean magnitude 
variations are indeed periodic. 

We estimated the relative scale of photometric period vari- 
ations within their 3cr limits to be Z = 0.020 for LQ Hya. 
Assuming that these variations are caused by starspots on dif- 
ferent stellar latitudes experiencing differential rotation, we de- 
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rive an estimate k ^ 0.02 for the differential rotation coefficient 
or equivalently AQ x 0.078 rad/d for the surface shear. This 
value indicates an almost rigid rota tor and is v irtually the same 
as the estimates bv lJetsul (Il993h and lYoul (l2007l) . Alternatively it 
is possible that the period variations are caused by factors unre- 
lated to surface differential rotation such as active region growth 
and decay or migration governed by the underlying dynamo field 
(Korhonen & Elstner 2011) and that the signal from any weak 
differential rotation gets buried under these. In either case our re- 
sult conforms to theoretical results whic h indicate that fast rotat- 
ing s tars should approach rigid rotators (iKitchatinov & Riidigeij 
[19991) . 

On the basis of the light curve minima, which indicate the 
main spot regions on the star, the typical configuration of the 
spot activity on LQ Hya seems to be one or two active regions at 
different longitudes. These active regions seem to be relatively 
stable in time scales up to six months, although they usually un- 
dergo longitudinal migration during that time. Quite often the 
two active regions have moved so close to each other that we 
only observe one merged light curve minimum. Thus, the ex- 
istence of even a third major active region producing a merged 
light curve minimum with either of the other two regions re- 
mains possible. 

Contrary to the relatively short term stability of the individ- 
ual active regions, there have been remarkably little stable struc- 
tures in the long term phase distribution of the light curve min- 
ima. We found one active longitude, with a possible secondary 
active longitude at Acp = 0.5 phase separation from the primary 
rotating with the period P^i = 1.603693 ± 0.000058 d. However, 
this active longitude has only been present in the data between 
the years 2003 and 2009. Before this we found one possible 
structure during 1995 and some time before that rotating with 
the period P^i = 1.68929 ± 0.00008 d. Because this period is so 
long compared to the other found light curve and active longi- 
tude periods, it is entirely possible that it does not correspond 
to any physical structur e on the star This picture is contrary to 
iBerdvugina et al.l (l2002l) . who report two persistent active longi- 
tudes for the total 20 years of their analysed photometry. We find 
no persistent active longitudes, nor any regular flip-flop events. 

We note that Pai = 1 .603693 d is close to the mean photo- 
metric period P„ - 1 .6043 d and well within the fluctuation of 
the individual light curve period estimates. In another active star 
II Peg there were indicat ions that the active longitude rotated 
faster than the star itself (iHackman et al.ll2012t iLindborg et al.l 
1201 ih . which could be explained by a dynamo wave propagating 
in the azimuthal direction (e.g. Krause & Radler 1980). There is 
no indication of such a dynamo wave in LQ Hya. 

The mean timescale of change for the light curve of LQ Hya 
is Tc = 50.5 d. This is longer than the maximum length of the 
individual datasets ATaax which means that the light curve typ- 
ically stays stable within the individual datasets. As a compari- 
son, we estimated the c onvective turnover time to be t^ - 33.5 
d using the formula of lOssendriiverl ("1997) (our Eq. |4). These 
two values are of the same order of magnitude. Previously we 
estimated a very similar ratio for these values, i.e. Tq = 44.1 
d and Tc = 28.5 d, f or the young solar analogue HD 116956 
(iLehtinen et al.ll201l1) . Although it is not clear whether there is 
an actual connection between the two values or if their ratio is 
simply governed by the choices made in the numerical proce- 
dure, such a result could be expected. After all, the strength of 
convection in the star is bound to be a strong factor in the evolu- 
tion of the spot structures on late-type stars. 
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